Astrophysics and Space Science manuscript No. 

(will be inserted by the editor) 



V.S.Beskin • E.E.Nokhrina 

On the role of the current loss in radio pulsar evolution 



o 
o 

(N 
CO 



CN 
> 
OS 

oo 

00 

o 
o 

6 

o3 



Received: date / Accepted: date 

Abstract The aim of this article is to draw attention to 
the importance of the electric current loss in the energy 
output of radio pulsars. We remind that even the losses 
attributed to the magneto-dipole radiation of a pulsar in 
vacuum can be written as a result of an Ampere force ac- 
tion of t he electric currens flowing over the neutron sta r 
surface (|Michell99lHBeskin, Gurevich fc Istominlll993ft . 
It is this force that is responsible for the transfer of angu- 
lar momentum of a neutron star to an outgoing magneto- 
dipole wave. If a pulsar is surrounded by plasma, and 
there is no longitudi nal current in its magnetosphere, 
there is no energy loss (IBeskin. Gurevich fc Istomin 1993; 
iMestel. Panagifc Shibatal Il999h 7 It is the longitudinal 
current closing within the pulsar polar cap that exerts 
the retardation torque acting on the neutron star. This 
torque can be determined if the structure of longitudi- 
nal current is known. Here we remind of the solution 
by Beskin, Gurevitch & Istomin (1993) and discuss the 
validity of such an assumption. The behavior of the re- 
cently observed "part-time job" pulsar B1931+24 can 
be naturally explained within the model of current loss 
while the magneto-dipole model faces difficulties. 
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1 Magneto-dipole loss 

The first idea to explain the energy loss of rad io pul- 
sars w as to consider the magneto-dipole radiation (jPacinil 
Il967[ ). Indeed, the magneto-dipole formula gives for the 
radiation power 
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where \ is t ne angle between rotational and magnetic 
axis, R is a neutron star radius ~ 10 km, and Q is a pul- 
sar angular velocity. This formula explains pulsar activ- 
ity and observed energy loss for expected large magnetic 
field near the surface Bo ~ 10 12 Gs. 

Let us recall that the physical reason of such en- 
ergy loss is the action of the torque exerted on the pul- 
sar by the Ampere force of th e electric curr ens flowing 
over the neutron star surface (|lstominll2005ft . The elec- 
tric and magnetic fields in the outgoing magneto-dipole 
wave in vacuum can be found by solving the wave equa- 
tions V 2 B + i? 2 /c 2 B = and V 2 E + n 2 /c 2 E = with 
the boundary conditions stated as the fields correspond- 
ing components Et and B n being continuous through 
the neutron star surface. Inside the star one can con- 
sider magnetic field as homogeneous, and find the corre- 
sponding electric field using the frozen-in condition. As 
a result, the full solution will give us the discontinuity of 
{B t } and {E n } attributed to the surface charge a s and 

the surface current J s = — [n. {B}1. The Ampere force 

47T 

exerts the torque 

K = ±J[r,[J„B]]dS (2) 

on the neutron star. The energy loss of a pulsar due 
to this torque is equal to Jl}. Thus, it is the surface 
current that is responsible for the angular momentum 
transform fro m a neutron s tar to an outgoing magneto- 
dipol e wave ( Michell 1 1 9 9 it : IBeskin. Gurevich fc Istominl 
11993ft . 

A pulsar in vacuum loses its rotational energy due 
to angular momentum transform to the electromagnetic 



2 



V.S.Beskin, E.E.Nokhrina 



wave at the rate given by ([I]). However, this is not so 
if the pulsar magnetosphere is filled with plasma and 
there is no long itudinal current in the magne t ospher e. As 
was shown by [Beskin. Gurevich fc Istomir] ( 19931 ) and 
iMestel. Panagi fc Shibatah 1999t ). in this case the Poynt- 
ing flux through the light cylinder is equal to zero. In- 
deed, as the ideal conductivity condition is applicable not 
only inside the neutron star but outside as well there is 
no magnetic field discontinuity at the star surface. Con- 
sequently, there is no Ampere force acting on a pulsar 
and, hence, there is no energy loss. For zero longitudinal 
current the light cylinder is a natural boundary of the 
pulsar magnetosphere. 



condition that no longitudinal current can flow in the 
region of the closed field lines. 

For arbitrary inclination angle x the current io can be 
written as a sum of its symmetric ig and anti-symmetric 
ix components. The anti-symmetric current begins play- 
ing the main role when the pulsar polar cap crosses 
the surface where the Goldreich-Julian charge density 
Pgj = O • B/2-7TC changes sign. This condition can be 
written as 



X 



7T OR 



(6) 



2 v c 

For example, taking the Goldreich-Julian current density 



2 Current loss 

In this section we remind the exact solution for the sur- 
face curr ent within the polar cap presented in the mono- 
graph bv lBeskin. Gurevich fc Istominl ( 19931 ). As it was 
shown in the previous section, the neutron star retar- 
dation is due to Ampere force Fa = J s x B/c. If the 
magnetosphere is filled with plasma, the surface current 
J s is flowing within magnetic polar cap only. This sur- 
face current closes the volume longitudinal current in the 
magnetosphere and the return current flowing along the 
separatrix between open and closed field lines region. 

In order to write the equation for the surface current, 
the several assumptions must be made. We assume that 
the conductivity of the pulsar surface is uniform, and 
the electric field E s has a potential, so that the surface 
current can be written as J s = V£'. Using the stationary 
continuity equation divJ = 0, where dJ z /dz is equal 
to the volume current %uBq flowing along the open field 
lines, one can obtain 



V 2 £' = -i N B 



(3) 



Making the substitution x = sm9 m and introduc- 
ing the non-dimensional potential £ = Att^' /BqR 2 S1 and 
current iq = — Airiu/QR 2 we get 
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Here 9 m and ip m are polar and azimuth angles with re- 
spect to magnetic axis. 

Equation ([4]) needs a boundary condition. This bound- 
ary condition results from the proposition that there is 
no surface current outside the magnetic polar cap. This 
means that 



£ [xo(<Pm), <Pm] = COnst, 



(5) 



where xo(ip m ) is the polar cap boundary. The solution 
for the Dirichlet problem fU [5|) can be obtained by the 
Fourier method. The jump in the potential derivative 
at x = xo(ip m ) gives us the current flowing along the 
separatrix. As we see, it is defined uniquely by the lon- 
gitudinal current in the region of open field lines and by 



iGj(x,tp m ) « cosx+-xcos(,9 m sinx = is+i-AX cos p m , (7) 

we obtain the following solutions of the Dirichlet prob- 
lem (]4])-([5]) for the symmetric and anti-symmetric volume 
currents respectively: 
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£a = ~^-x(x - X ) COS <p r , 



(8) 
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The torque exerted by the surface current over the 
neuron star can be written as 



K = - 

c 



- c J [r, [J„(B )]]dS, 



(10) 



where Bo is the dipolc field near the neutron star surface. 
Let us decompose the braking torque K over the orthog- 
onal system of unit vectors e m , ni, and ri2. Here e m is 
a unit vector along the magnetic moment; vector ni is 
perpendicular to the magnetic moment and lies in the 
plane of the magnetic moment and the rotational axis; 
vector ri2 complements these to the right-hand triple: 



K = Kue m + K ± n 1 + K f n 2 . 



(11) 



i (x,ip m ). (4) a 



iff plays no role in Eulcr equations that describe the 
rotational dynami cs of the decelerating neutron star . As 
result we have (jBeskin. Gurevich fc Istominl 1 19931 ) : 
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where A = x cos ip m d^/dx — sincp m d£/d(p m . 

The leading perpendicular torque component K\ is 
equal to zero equivalently for arbitrary shape of the polar 
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cap due to the boundary condition ((5|). The values K\ 
and Kj_ can be written as 

1/2 



K \\ = 
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(16) 
.(17) 



Here the coefficients and /i^ depending on the shape 
of the polar cap are much less than unity, and the coef- 
ficients ci and c±_ ~ 1. 

We can now find the derivatives of the angular ve- 
locity fl and of the inclination angle x of a neutron star 
through the Euler dynamics equations: 



if_Lsinx, 



(18) 
(19) 
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~dt = 11 C ° SX 
J r f2—^- = K± cosx — K\\ sinx- 

For the inclination angles x n °t too close to 90° (i.e., for 
cosx > (^-R/c) 1 ^ 2 ); when the anti-symmetric current 
plays no role in the neutron star dynamics, we obtain 



dt 

&X 
dt 
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is cosx, 



B%f2 2 R 6 
C|| i s sin x- 
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-z s cosx- 



As a result, for homogeneous current density within 
open magnetic field lines region is = ju/jGJ = const 
where jqj — cpQj we have 

4 c 

Here /* is the non-dimensional area of a pulsar polar 
cap: S'cap = f*Tt(f2R/c). It depends on the structure of 
the magnetic field near the light cylinder. For a pure 
dipole magnetic field (and aligned rotator) /* = 1, and 
for a magnetosphere containing no longitudinal currents 
changes fr om 1.592 for the aligned rotator, x = 0° 
(iMichelj 199lh. to 1.96 for an orthogo nal rotator, x = 
90° (|Beskin. Gurevich fc Istomin| [T993). Recent numeri- 
cal calculations for an axisymmetric magnetosphere with 
non- zero longit udinal electric current give f* sa 1.23 - 
1.27 (|Gruzinovll2005l : lKomissarovll2006l : lThnokhinll2006h . 
If the singular point separating open and close field lines 
can be located inside the light cylinder, the value /* can 
be even ^> 1. As the Goldreich- Julian charge density near 
the polar cap is propotrional to cosx, one can write 

„ 7 /* Bq&r 6 2 

4 c J 
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3 PSR B1931+24 



The recent discovery of the " par t-time job" pulsars PSR 
B1931+24 (jKramer et al.l l2006l ) with f2 n/^ s ~ 1.5 
shows that the current loss is indeed playing an impor- 
tant role in the pulsar energy loss. If we assume that in 
the on-state the energy loss is connected with the current 
loss only and in the off-state with the magneto-dipole ra- 
diation (in which case the magnetosphere must be not 
filled with plasma) we get 



3/* 2 
~Y cot 2 x. 



(25) 



It gives x ~ 60°. On the other hand, if we assume 
the Spitkoysky's relation for the on-state energy loss 
flSDitkovskj2006h 
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we obtain 



f2 on 3(1 + sin x) 
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(26) 



(27) 



Clear, this ratio cannot be equal to 1.5 for any incli- 
nation angle. This discrepancy can be connected with 
that fact that all the numerical calculations produced 
recently contain no restriction on the longitudinal elec- 
tric current magnitude. As a result, current density can 
be much larger than Goldreich-Julian one. 



(22) 4 On the magnitude of a surface current 



As we have shown, the current loss plays the major role 
in the pulsar dynamics. In particular, the behaviour of 
the pulsar B1931+24 can be naturally explained within 
this model. The current loss model have some important 
consequences: 

1. the energy loss of an orthogonal rotator is QR/c 
times smaller than for an aligned rotator. This is con- 
nected with the boundary condition ([5]) that leads to 
almost full screening of the tor oidal magnetic field 
i n the open field lines region (see iBeskin &; Nokhrinal 
(Hoi)); 

2. consequently, during its evolution a pulsar inclination 
angle tends to 7r/2 where energy loss is minimal. 



On the other hand for x ~ 90° when the anti-symmetric 
current plays the leading role we obtain 



Wr. 
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(24) 



As we see, the energy loss of the orthogonal rotator are 
flRfc times smaller than of the aligned rotator. 



On the other hand, it is known for the Michel's monopole 
solution that in order to have the MHD flow up to in- 
finity, the toroidal magnetic field must be of the same 
order as the poloidal electric field on the light cylin- 
der. If the longitudinal current jii does not exceed by 
{QR/c) -1 / 2 times the respective Goldreich-Julian cur- 
rent density (for the typical pulsars this factor approach 
the value of 10 2 ), the light surface |E| = |B| for the or- 
thogonal rotator must be located in the vicinity of the 
light cylinder. In this case the effective energy conversion 
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5 Conclusions 

As we have seen, the current loss connected with the lon- 
gitudinal current flowing in the magnetosphere plays the 
main role in pulsar dynamics, and recent observations of 
"part-time job" pulsar supports this point. This evolu- 
tion includes not only a neutron star retardation but also 
the sufficient change in the angle between magnetic and 
spin axis. We have seen as well that the model of current 
loss depends crucially on the distribution of the electric 
current and its value in the inner gap. For current loss 
model with j'gj the inclination angle grows with time so 
a pulsar tends to be an orthogonal rotator. In this case 
the energy loss is to be fiR/c times smaller than for the 
aligned rotator. As a consequence, the light surface must 
be located in the very vicinity of the light cylinder. On 
the other hand, to realize the homogeneous MHD out- 
flow up to infinity for the orthogonal rotator the current 



